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Abstract 



According to resummed perturbation theory, certain scalar theories have a global sym- 
metry, which is restored in the vacuum but is broken at high temperatures. Recently, 
this phenomenon has been studied with 4d finite temperature lattice simulations, and 
it has been suggested that the non-perturbative dynamics thus incorporated would 
hinder the transition. We have carried out another lattice study, for a theory with very 
small coupling constants. We find perfect compatibility with next-to-leading order re- 
summed perturbation theory, and demonstrate that "inverse" symmetry breaking can 
indeed take place at high temperatures. 
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1 Introduction 



Inverse symmetry breaking is a phenomenon where a (global) symmetry is restored in 
the vacuum, but is broken at high temperatures. This is in contrast to the usual be- 
haviour, where symmetries may be broken at low temperatures but get restored at high 
ones. Originally, inverse symmetry breaking was shown to be possible in 0(N)xO(M) 
scalar theories 0, but something analogous might happen in more realistic cases as 
well. Indeed, several potential cosmological consequences have been discussed (see, 
e.g., 0] and references therein). 

Recently, there has been some interest in studying inverse symmetry breaking with 
lattice simulations || |]]. This interest is based on controversial statements existing 
in the literature, concerning the possibility of inverse symmetry breaking when non- 
perturbative effects are taken into account (see, e.g., 0-0 and references therein). The 
result of the simulations, both in three dimensions (3d) [[3] and four dimensions (4d) M , 
was that no sign of inverse symmetry breaking was seen, even though perturbatively 
it was supposed to take place. 

The purpose of this paper is to study inverse symmetry breaking with 4d finite tem- 
perature lattice simulations, in a theory that is coupled weakly enough for perturbation 
theory to work. We do observe clear signs of inverse symmetry breaking. We also dis- 
cuss the general way in which non-perturbative effects may manifest themselves in 
scalar theories of the studied type, by constructing the relevant dimensionally reduced 
3d effective theory. 

To be more specific, consider the Euclidian scalar Lagrangian 

£cont = i(9 M 1 ) T (9 M 1 ) + i(9 M 2 ) 2 + -m^! + -m^ 

+ ^Ax(^) a + ^A 2 ^ + i 7 ^i^, (1) 

where <pi is an O(N) symmetric real vector (we take N= 4) and 2 is a Z 2 symmetric real 
scalar field. The symmetry groups appearing here are not essential for the existence 
of inverse symmetry breaking. In |3|, |J, a case where both fields are real scalars and 
the symmetry is Z 2 X Z 2 was considered, but we have chosen to take an 0(4) field to 
undergo inverse symmetry breaking, since the finite volume effects in the simulations 
are then somewhat easier to control (a ^-field in the broken phase has problematic 
"tunnelling correlations" , because of a barrier separating the two degenerate minima) . 

Consider now the allowed values of the coupling constant 7 in eq. ([!]). In principle, 7 
can be negative, even though not by an arbitrarily large amount. This is because the 
stability of the potential at zero temperature requires that Ai 0, A 2 ^> 0, AiA 2 3> 97 2 . 
However, within these requirements, it may happen that (N + 2)Ai + 37 <C or 
A 2 + A^7 <C 0. These combinations of the coupling constants turn out to multiply 
the temperature squared in the effective finite temperature mass parameters of <pi and 
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02, respectively. Thus, at high temperatures, one of the symmetries (but not both of 
them) may get broken, even if both symmetries are restored at zero temperature. This 
is called inverse symmetry breaking. If a symmetry is broken already at T = 0, one 
talks of "symmetry non- restoration" . 

The obvious way to establish the phenomenon of inverse symmetry breaking would 
be to fix m\ > so that <pi is in the symmetric phase at T = 0, and then to increase 
the temperature. Equivalently, one can fix T and tune m\: it is then sufficient to show 
that there is a phase transition at some critical value m\ c > 0. Although we will use 
4d finite temperature simulations to study inverse symmetry breaking, we will use a 3d 
effective field theory to predict the values of the parameters where this phenomenon 
should occur. 

In the next section, we derive the perturbative estimates for inverse symmetry break- 
ing in the scalar theory of eq. (|1|) in some more detail, using the method of dimensional 
reduction and 3d effective field theories ||. This allows us to implement the resum- 
mations needed at finite temperature in a transparent way and, furthermore, to see 
what kind of non-perturbative effects there can be. In sect. 3, we study the theory in 
eq. ([TJ) with 4d finite temperature lattice simulations. The results of the simulations, a 
comparison with perturbation theory, and our conclusions are in sect. 4. 



2 The 3d effective theory: 

perturbative and non-perturbative results 

It is well known that perturbation theory at finite temperature requires resumma- 
tions. A convenient way to implement the resummations is to construct an effective 
3d field theory for the zero Matsubara modes. The construction of the 3d theory is 
purely perturbative, as only massive degrees of freedom are integrated out. The non- 
perturbative dynamics of the system is contained in the final effective 3d theory; it can 
thus be studied in an economic way with this approach. 

In more concrete terms, the purpose of this section is to derive the coefficients ci, c 2 
in the expression 

^ = Cl |7l+c 2 |7| 3/2 + C3|7r + ..., (2) 

c 

where T c is the temperature where inverse symmetry breaking takes place, and we 
assumed that parametrically Ai,A2 ~ and m\ ~ 0. Moreover, we review why, in 
this expression, non-perturbative effects only affect the coefficient C3 and the higher- 
order terms. In general, non-perturbative effects also affect the order of the phase 
transition (the transition is of second order, as we will see), but the existence of the 
transition itself (i.e. the coefficients 01,02) can be deduced purely perturbatively. 

The general method of dimensional reduction, its accuracy, and generic rules for 
its application M, have been discussed extensively in the literature (see also [HI ITT 
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and Jl2| for a review). Hence we discuss only the results here. 



In order to derive the coefficients ci, c 2 , it is sufficient to perform dimensional reduc- 
tion at the tree-level for the coupling constants and at the 1-loop level for the mass 
parameters. In the body of the text we work at this order, but in the appendix we 
give also the next corrections (of the relative order 0(|7|) with respect to the leading 
terms). This serves to show how the MS scheme scale parameter appearing in the 
zero temperature parameters gets fixed. Moreover, the 2-loop corrections in the mass 
parameters are needed if the non-perturbative constant C3 were to be determined with 
lattice simulations in the 3d effective theory. 

The first step of dimensional reduction is the integration out of non-zero Matsubara 
modes. After rescaling the fields squared and coupling constants by T, the final effective 
3d theory is 



S cS = / d 3 *^ WifWi) + -(«9,0 2 ) 2 + -m? >3 #>i + -m\^\ 

+ ^A li3 (0f0l) 2 + ^2^+^730^!^}. (3) 



To this order, the fields appearing are just the same as the zero Matsubara components 
of the 4d fields, apart from the trivial rescaling with T. To leading meaningful order, 
the expressions for the parameters appearing in eq. (^) are 

Ai, 3 = TAi, A 2i3 = TA 2 , 7 3 = T7, (4) 

< 3 = m ? + ^( v ^— A i+7j, ml 3 = ml + — \\ 2 + N>y J, (5) 

where the parameters on the RHS mean renormalized MS scheme parameters at some 
scale ~ T (to be more precise, see the appendix). 

Near the symmetry breaking phase transition, the theory in eq. (§) can be further 
simplified. Indeed, recall that we have chosen it to be the O(N) field 0i that undergoes 
inverse symmetry breaking: 

(A^ + 2)A! + 37<0. (6) 

Around the critical temperature, m 2 3 ~ (more precisely, as we will see, m 2 3 ~ 
|7| 3 / 2 T 2 ). On the other hand, it follows from eq. @ and the vacuum stability require- 
ment A1A2 ^> 97 2 that 

A 2 + A^ 7 > 0. (7) 

Hence the effective mass parameter m\ 3 is "large", m\ 3 ~ \l\T 2 ^> wf 3 , and the 
heavy excitations corresponding to the field 2 can be integrated out. 
The final effective 3d theory after integrating out 2 is of the form 

S eS = /^{^(9A) T (^0i) + iml 3 0r0i+ 2 \Ai,3(0f0i) 2 }. (8) 
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The parameters that appear are 



Ai,a = A li3 , (9) 

/ _/y _|_ 2 \ 1 

™i, 3 = m l + 24 l,^ - Al + 7 J ~ 8^ 7sm2 ' 3 ' ^ 10 ^ 

The next corrections are given in the appendix. 

Using these expressions, we can discuss the perturbative (and non-perturbative) 
predictions for inverse symmetry breaking. To be specific, let us fix the parameter 
values. To safely satisfy eq. (f|) and the vacuum stability requirement, we choose 

T = -iTl, Ai = 2(jV 3 +2) l7|, A 2 = 10(iV + 2)| 7 |. (11) 

These refer to the MS scheme parameters at some scale /io- Note that even though 
they are parametrically of the same order of magnitude, the coupling constants in 
eq. (|TTD are numerically widely different. In particular, m\ z involving X 2 is "large", 
and therefore the first correction it induces in m\ 3 (the last term in eq. ( |IDD ) is also 
significant. However, further corrections are relatively much smaller, by (73/^2,3). 
To further simplify matters, let us take m\ ~ 0. 

The phase structure and critical temperatures. We have derived above the 
effective 3d theory in eq. (§), and argued that it describes the infrared properties of 
the system non-perturbatively around the point m\ 3 ~ 0. What does this imply? 
The non-perturbative properties of the theory in eq. @ are of course well known: the 
theory has a second order phase transition into a broken phase when m\ 3 becomes 
negative. The exact value of the MS scheme parameter m\ 3 at the phase transition 
point could be determined with lattice simulations, but this is not important for us 
here. It is sufficient to know that, since the coupling Ai,3 is dimensionful, the result 
can only be of the form 77^ 3^ = (c/lQir 2 )Xf 3 , where c ~ 1 (c could also be negative). 
The factor l/(167r 2 ) arises from typical 2-loop order (see the appendix) contributions 
to c, cancelling the scale dependence of m\ z . 

Using now the expression in eq. (0); the value of m 2 ^ 3 rit ' ) can be converted to the 
critical temperature. It is seen that the non-perturbative constant c only contributes 
to c 3 in eq. (^). 

Thus, the constants Ci, c 2 in eq. (|2|) can be derived simply by equating eq. ( |TU| ) with 
zero. This corresponds to the next-to-leading order accuracy for critical temperature 
discussed in general in |l3j , and in the present context in || . Using eq. (fLl]^ ,T ' - : 
becomes 



m I,3 



_n\_ |7 



3 / 2 20 + 9N 



:i2) 



48 8tt V 24 

To have symmetry breaking, the coefficient of T 2 must be negative, so that we need, 
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for iV = 4, t 

24 




, r = 0.34. (13) 
20 + 9iV v ' 

To satisfy this requirement, we fix in practice |7| 1//2 = 1/6. Then the critical tempera- 
ture for inverse symmetry breaking is determined by 




ml_ H H 3 ' 2 ./20 + 9iV K3xl0 ^ 



T} 48 8tt V 24 



The scalar field expectation value. Besides the critical temperature itself, we 
would also like to know the behaviour of around T c . It turns out that a good 
enough estimate can already be obtained with the tree-level result in the effective 
theory of eq. (W): 



(<M = \ my,*,. A ; „ (15) 

We have also computed the 2-loop effective potential in the theory of eq. @, as can 
be easily done (then one also has to use the 2-loop parameters in the appendix), but 
numerically the effects found are very small for the small coupling constants used, 
typically ~ 1%. 



3 Finite temperature 4d lattice simulations 

We now turn to 4d finite temperature lattice simulations, and attempt to verify 
eqs. (0), (0)- As a first point, let us discuss the renormalization of the theory. 

As mentioned above, the coupling constants considered are very small. Thus pertur- 
bation theory sensitive only to ultraviolet (UV) degrees of freedom is well convergent. 
As a consequence, the theory can be renormalized in perturbation theory. 

To be more specific about the renormalization, recall that in sect. |]we have param- 
etrized the theory by couplings defined in the MS scheme at an arbitrary scale /2 . 
In the simulations, in contrast, one uses the lattice regularization. To compute the 
relation between the schemes is a standard exercise (it is analogous to computing the 
relation between Aqcd in the continuum and on the lattice ||14|| , but much simpler as 
this is a scalar theory). One computes a set of 1-loop graphs and requires that the 
results are the same in the MS and lattice regularization schemes. At the 1-loop level, 
the result is as follows. Let 

L a (p) = — — \2 ln(o/2) + 167r 2 ri + ll . (16) 



Then the bare parameters appearing in the lattice action are 

2 2 (iy + 2)A 1 + 3 7 r (iV + 2)Aim[+37m| 
m 1B = m^fio) L a (fi Q ), 
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2 \ 2 + N<yr A 2 m 2 , + iV 7 m 2 

m 2B = «Wo) ~ 2 ~ L a (no), 

a z 2 



AiB = Ai(/i )-- — - — A x + 7 L a (n ) 



2 V 9 

A 2B = A 2 (/x )- 2 (A2 + iV 7 2 )L a (/io), 

7b = 7(^o) - ^ (4 7 2 + iV ^" 2 7Ai + 7 A 2 ) L a (A))- ( 17 ) 

These bare parameters are of course independent of JIq. The constants that appear 
are 

r = 0.154 933 390, n = -0.030 345 755. (18) 
The lattice spacing appearing here is determined by 

' (19) 



N+T' 



where N t is the temporal extent of the lattice. 

When scaled into a dimensionless form by fa — > (y/Kl/ a) fa, the lattice action becomes 



S = -Ki ^24>i(x + fi)fa{x) - K 2 ^2fa(x + fi)fa{x) + ^2<f>[(x)fa(x) + Y^^l( x ) 

X,fJ, X,fl X X 

+ ^[^^W-lf + ^^^-lf + a^^M^^W. (20) 

X XX 

Here the parameters are related to those in eq. ([17]) by 
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K>i — — 
Kb 



-(4 + \a 2 m 2 B ) + J(4 + ^%? B ) 2 + ^A iB 



A = AiB«? , a = ^7s«i«2- (21) 

The general philosophy of the simulations is now as follows. It is seen from eq. ( |T2"D 
that, for a fixed m 2 , the O(N) symmetry is expected to be broken when T is increased 
above some critical value. Equivalently, as mentioned, one can fix T and tune m\ : when 
m\ is smaller than some critical value, but still positive, there should be symmetry 
breaking. The latter viewpoint is easier to implement in simulations since, according 
to eq. (|I9"D, the lattice spacing is then kept fixed. As there is only one temperature, we 
can also choose, for convenience, the arbitrary fixed parameter Jxq to coincide with T. 
Then, we just vary m 2 (/i )/T 2 . 

The parameter values following from eqs. fllTf) , with coupling constants as in eq. fllT ) 
and mKfio) = 0, are shown in Table [I] together with the lattice sizes used. The 
numerical simulations have been performed at small values of f3\ and f3 2 in eq. (|20"D 
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mj/T 2 Ki/IQ- 1 K2/IO- 1 /3i/10- 5 /3 2 /10" 3 a/10" 4 lattices 



5.0xl0" 5 

1.5xl0~ 4 
3.0xlO" 4 
6.0xl(T 4 
1.0xl0~ 3 
2.0xl0~ 3 
3.0xl0- 3 
4.0xl0~ 3 
6.0xl0~ 3 



1.2497845 

1.2497806 
1.2497747 
1.2497630 
1.2497474 
1.2497083 
1.2496693 
1.2496302 
1.2495522 



1.2570608 

1.2570608 
1.2570608 
1.2570609 
1.2570609 
1.2570610 
1.2570610 
1.2570611 
1.2570613 



N t = 2 
1.8185293 

1.8185179 
1.8185009 
1.8184668 
1.8184213 
1.8183077 
1.8181941 
1.8180804 
1.8178532 
N t = A 



4.7496966 

4.7496967 
4.7496968 
4.7496969 
4.7496972 
4.7496977 
4.7496982 
4.7496988 
4.7496999 



-4.4763402 

-4.4763263 
-4.4763053 
-4.4762634 
-4.4762076 
-4.4760680 
-4.4759283 
-4.4757887 
-4.4755096 



{12 3 ...36 3 } 

12 4 ...20 4 

{20 3 ,24 3 } 

{20 3 ,24 3 } 

{20 3 ,24 3 } 

{20 3 ...28 3 } 

{20 3 ,24 3 } 

{20 3 ...28 3 } 

{20 3 ,24 3 } 

{12 3 ...36 3 } 



5.0xl0~ 5 1.2497859 1.2568207 1.8213989 4.8442880 -4.5055481 {24 3 } 
6.0xl0~ 3 1.2497278 1.2568209 1.8212295 4.8442891 -4.5053391 {24 3 } 



Table 1: The coupling constants and the lattice sizes used. The numbers in the 
curly brackets indicate the spatial volume; the temporal extent is N t . The notation 
Z/™ in ...L™ ax means that the lattice extent L is increased in steps of 4 from L min to L max . 



where the cluster algorithm is not very effective. We have chosen instead a hybrid over- 
relaxation algorithm |16[]. For each point we performed 50K sweeps for thermalization 
and 100K sweeps for measurements. The errors were always computed from a blocking 
analysis searching for a plateau behaviour. 



4 Results and conclusions 

In Fig. [TJ, we show a scan of continuum parameters, for different volumes. The field 
is shown in lattice units, 1 1 att = (a/ ^Ki)<f)l ont . The infinite volume curve computed 
from the 2-loop effective potential in the theory of eq. (§) is also shown (however, as 
mentioned, the tree-level result could also be used, as the difference is very small, of 
order 1%). 

It can be seen that the lattice results differ from the continuum perturbative re- 
sult by quite a significant amount for the volumes used. However, this can be easily 
understood. Indeed, the finite volume mean field estimates obtained from 



/c^expt-liVWeff^i)] 
JdMleM-T N3a3V es(4>i)} 



(22) 



where V e g((j)i) is the continuum effective potential, 0i denotes the radial component 
of 0i, and the volume element of an 0(4) vector is d<f>\ oc 3 <i0i, are also shown in the 
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Figure 1: A scan of lattice results, as a function of the continuum parameter m 2 /T 2 . 
The infinite volume perturbative continuum result, as well as finite volume mean field 
estimates, are also shown. The lattice results are completely compatible with the 
perturbative ones. 



figure. They are in perfect agreement with the lattice results. 

Next we take two points and look more precisely at the approach to the infinite 
volume limit. The results are shown in Fig. 0. Again we find perfect compatibility 
with the perturbative mean field estimates. For the parameter value mf/T 2 = 0.00005, 
there is inverse symmetry breaking, while for m 2 /T 2 = 0.006, the symmetry is restored 
in the infinite volume limit. 

The results so far had been obtained with a single lattice spacing, N t = 2. We have 
also made simulations with N t = 4. The results are shown in Fig. [| after rescaling 
volumes and units so that they become comparable with those at Nt — 2 (in other 
words, the volume 4 x 24 3 is represented at the same point as 2 x 12 3 , and the field 
has been normalized so that the data points correspond to the same continuum units). 
No lattice spacing dependence is seen within the errorbars. 

Finally, let us point out that, in principle, the existence of inverse symmetry breaking 
could also be verified without any reference to perturbative renormalization. In prac- 
tice, though, this would require huge volumes for the weak couplings we have chosen. 
However, we have made a few zero temperature runs for the parameters corresponding 
to m\/T 2 = 0.00005. The results are shown in Fig. ||] with the filled circles. It is seen 
that the results for are indeed smaller than at finite temperature, and certainly 
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Figure 2: The infinite volume limit for the 0(4) field expectation value (<fi\) at m\/T 2 = 
0.00005, 0.006. Mean field estimates are also shown, as well as results from a smaller 
lattice spacing. 



consistent with an extrapolation to zero. 

In conclusion, we have reiterated what kind of non-perturbative effects there can 
be in finite temperature phase transitions in scalar theories. The existence of the 
transition alone can be deduced with perturbation theory, while some of its properties 
(such as the order) are non-perturbative. With explicit 4d finite temperature lattice 
simulations, we have verified that the non-perturbative behaviour is indeed in perfect 
agreement with perturbative estimates, and inverse symmetry breaking takes place. 
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Appendix 

We give here the parameters appearing in the 3d theories in eqs. (^) and (^) with more 
precision than in the text. 
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Denoting c = -0.348725 and 



L b {p) = 2 In 



2 In A 



47re-TT /i T 
the MS scheme renormalized parameters appearing in eq. (|3|) are 

3 /iV + 8 



(23) 



Al,3 
^2,3 
73 



= T 
= T 
T 



AiO) 



-a 2 + 7 2 )^) 



32vr 2 V 9 
X 2 (fl)-J-(xl + N 1 2 )L b {fl) 



70) 



32vr 2 
1 



N + 2 



, v 4 7 2 + ^^7Ai + 7A2) L b (p) 
T f N + 2 



32tt 

mi(/i) -32^l^- Aimi+7m2 J L6(/i) + 24V 3 



Al,3 + 73 



1 (N + 2 



32tt 5 



9 ^,3 + 7ij^ y +c 



3T 



m^(/2) - — ^ ^A 2 m 2 + iV 7 m 2 )L fe (/i) + ^ ( A,., A \, 
1 ^A 2 , 3 + iV 7 3 2 )(ln^ + C 



32tt 2 V3 



(24) 



The 3d coupling constants here are scale independent, whereas the mass parameters are 
not (due to the terms on the latter lines), corresponding to the 2-loop UV-divergences 
remaining in the super-renormalizable 3d theory. When computing, for instance, the 
2-loop effective potential in the 3d theory, this scale dependence gets cancelled. Where 
not indicated explicitly, the parameters appearing in the expressions in eqs. (|4]) are 
renormalized MS scheme quantities. 

The parameters appearing in the final effective theory in eq. (|8]) are 



Al,3 



^1,3 (#) = m l(fiT) + 



3 7! 



16-k m 2 ,3 
2,-. n , T /N + 2 



24 V 3 Al ' 3 + 73 ) ~ S^ 73 ™ 2 ' 3 



+ 



32tt 2 



^ 2 A 3T 1 

^2,373 - 73 ^^ + 0+- 



N + 2 
9 



3T 



A?3(ln^ + c 



(25) 



.(26) 



For mf 3 , the first line is the 1-loop expression and the second line the 2-loop contribu- 
tion. In these formulas, it is enough to use the 1-loop part of m^, given in eq. (||). 
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